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ABSTRACT 

The Milky Way appears to be missing baryons, as the observed mass in stars and gas is well below 
the cosmic mean. One possibility is that a substantial fraction of the Galaxy's baryons are embedded 
within an extended, million-degree hot halo, an idea supported indirectly by observations of warm gas 
clouds in the halo and gas-free dwarf spheroidal satellites. X-ray observations have established that 
hot gas does exist in our Galaxy beyond the local hot bubble; however, it may be distributed in a hot 
disk configuration. Moreover, recent investigations into the X-ray constraints have suggested that any 
Galactic corona must be insignificant. Here we re-examine the observational data, particularly in the 
X-ray and radio bands, in order to determine whether it is possible for a substantial fraction of the 
Galaxy's baryons to exist in ~ 10 6 K gas. In agreement with past studies, we find that a baryonically 
closed halo is clearly ruled out if one assumes that the hot corona is distributed with a cuspy NFW 
profile. However, if the hot corona of the galaxy is in an extended, low-density distribution with a 
large central core, as expected for an adiabatic gas in hydrostatic equilibrium, then it may contain up 
to 10 11 M Q of material, possibly accounting for all of the missing Galactic baryons. We briefly discuss 
some potential avenues for discriminating between a massive, extended hot halo and a local hot disk. 
Subject headings: Galaxy: halo — X-rays: diffuse background 



1. INTRODUCTION 



The question of whether or not significant reservoirs 
of hot baryons exist aroun d bright field galaxies remains 
a topic of current debate ([Dai et al. 2010| Anderson fe 
Bregm an 2010; Humphrey et al. 2011| |Prochaska et al. 
2011 ) despite years of discussion in the lii 
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literature ( Spitzer 



Spjtzer| [19691 |Mo fc Miralda-Escude 
Theoretical prejudice favors the 
idea that quasi-stable ~ 10 6 K coronae develop as sh ock- 
heated material in the aftermath of halo collapse (|Silk 
19771 |White fc Rees||i~978t |White fc Frenk||199ll |Benson 
eFaL||2000| |Keres et al.l|2005[ |Dekel fc Birnboim| |2006[ 
Keres fc Hernquist||2009] |Crain et al.||2010p, though the 



density structure and mass content of these halos is sen 
sitive to uncertain physics and energy injection processes 
'White & Rees||1978| Mailer & Bullock||2004[ iKaufmann 



et al ||2U09[ |Bmsonp0T0[ |Sharma et al ||2U12| ). 

it hot gaseous halos exist, they would provide a poten- 
tial hiding place for missing galactic baryons - the > 80% 
of baryons that are un accounted for by colla 
Fukugita et"aL||l_998[ 
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It is not known whether the missing halo baryons exist 
pri marily as diffuse hot gas halos around normal galax- 
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ies (|Mallei_fc Bullock ||2004 
Sommer-Larscn 2006[) 

result of en ergetic blow out (|Dekel fc Silk|1986| |Almeida 



jeen mostly expelled as a 



et al. 2008 1 , or were never accreted i n the first place, pos- 
sibly as a result of pre- heating (e.g. Mo et al.||2005|). So 



far, searches have failed to detect exten ded X-ray emis- 
sion around nearby spiral galaxies (e.g., Bregman et al. 
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1998||Benson et al.|[2000| [Rasmussen et al.||2009| ). Most 
of the detected X-ray emission in these galaxies is cen- 
tered around disk/bulge region, and is associated w ith 



active star formation regions (see, e.g., Wang 2007 Li 



et al.||2007[ 

Locally, the Milky Way provides an important bench- 
mark for understanding the missing galactic baryon 
problem. The Galaxy's dark matter halo mass is some- 
what uncertain, but maser observations, stellar halo trac- 
ers, and satellite kinematics suggest a t otal virial masa^ 
in the range M v = (1 — 2) x 10 12 M Q (see Boylan-Kolchin 
et al. 2012 for a summary), within an associated virial 
radius K v ~ (260 — 330) kpc. In the absence of mass loss, 
the Milky Way's baryonic allotment therefore should be 
Mb = fbM v ~ (1.65— 3.3) x 10 11 M , assumin g a univer- 
sal baryon f raction = 0.165 (consistent with Komatsu 



et al. 2011). The observed cold baryonic mass of the 



Milky Wa y is well be low this, M* ~ 0.65 x 10 11 McMil 



lan & Binney (2012), with effectively negligible co ntri- 
butions from coo l disk gas (|Kalberla fc Kerp||2009"| and 
satellite galaxies ( McConnachie|2012[ ). At least 10" M 
of baryons are missing from the Galactic census. 

There is circumstantial evidence that at least some 
of these missing Galactic baryons are in an extended, 
hot corona. For example, a confining hot mediu m can 



help explain gas clou ds in the Magellanic Stream ( Stan 
imirovic et al.||2002|), and the shape s of shells along the 



edge of the Large Magellanic Cloud ( de Boer et al.|1998 1 
are similarly suggestive. H- alpha emission at the lead- 
ing edge of the Magellanic stream is best understood via 
heating from a fairly dense hot medium (~ 10 -4 cm 



at ~ 50 kpc; |Weiner fc Williams|[l996l ). The thermal 



pressure of high- velocity clouds at the same distance is 

4 We define the virial mass using an average enclosed density of 
A = 95 times the critical density of the Universe, as determined 
using the spherical top-hat collapse model in the LCDM cosmology. 



2 



Fang, Bullock, & Boylan-Kolchin 



consi stent with a similar confining medium (Fox et al. 
20051. More recent ultraviolet and optical absorption 
probes of the Magellanic Stream provide strong evidence 
for a complex, multi-phase structure that is being evapo - 
rated by interactions with a hot corona ( Fox et al.|2010 ). 
The head-tail structure of high velocity clouds is another 
observa tion best explained in the context of a confining 



The third case explores a scenario where all of the X-ray 
gas detections are explained by a local distribution with 
globally negligible mass content (< 10 8 M Q ). 

2. MODEL DEFINITIONS 

In our fiducial explorations, we assume that the Milky 
Way's dark matter halo has a virial mass of M v = 
corona ( TPutman et al.||2011j ). A separate argument for 10^M© (with associated radius R v = 260 kpc) and that 



the existence of such a medium is the lack of detected it follows an NFW profile with a concentra tion of C v = 12 
HI in most dwar f spheroidal galaxies within 270 kpc of 
the Milky Way ( |Grcevich fc Putman||"2009[ ) - a result 
that is naturally explained by ram pressure strip ping by 



(and thus a scale radius R s = 21.7 kpc; Navarro et al 



10 4 cm 3 hot medium at ~ 70 kpc (see als o M oore 
fc Davis||1994| [Lin fc Faber|l983| |Nichols et aL|2(Jll| ) 



While these indirect probes are enlightening, the most 
direct avenue for detecting a hot corona is via X-ray stud- 
ies. It is well known that the soft X-ray background (be- 
tween 0.1 and 2 keV) consists of three components: the 
Local Hot Bubble (LHB), extragalactic emission (mostly 
from active galactic nuclei), and a thermal com ponent 



that lies at an in termediate distance (see, e.g., Kuntz 



& Snowden 2000). One popular interpretation of this 
thermal component is that it originates from the hot 
interstellar medium (ISM) in the Galactic disk, i.e., a 
hot gas disk rather than a ~ 100 kpc extended Galac- 
tic halo. Indeed, us i ng join t X-r ay emission-absorption 



analysis, Yao et al. (2009a) and Hagihara et al. (20101 



argued that the hot gas m our Galaxy is confined within 
a few kpc around the stellar disk. The detection of the 
highly ionized metal absorption lines at z — in numer- 
ous quasar spectra also suggested a hot gas component 
in and around ou r Galaxy (see, e.g., Nicastro et al.|200"2l 
Fang et al.|f2003l [20061 |Bregman fc Lloyd-i)aviesP007[ 



Gupta et al. 2012p ; however, it is also unclear this ab 



sorpti on is produced by hot gas in the disk or distant 



halo. Anderson & Bregman (2010) focused specifically 



on the question of baryonic closure, asking whether the 
dispersion measure of pulsars in the Large Magellanic 
Cloud (LMC) could be reconciled with a baryonically 
closed Milky Way halo. They concluded that it could 
not, though they focused mainly on cases where the hot 
gas follows a cuspy density profile characterized by either 
a single power law or a broken power law (as expected 
for dissipationless dark matter). 

In what follows, we re-examine the question of the the 
Milky Way's hot corona in the context of X-ray surface 
brightness constraints and pulsar dispersion measure- 
ments from the standpoint of three illustrative models 
for the hot gas density distribution: (1) an extended pro- 
file with a central core, as expected for an adiabat ic gas 
in hydrostatic equilibrium ( |Maller fc~ Bullock 2004 , here- 
after MB) ; (2) a centrally concentrated |jNavarro, Frenk, 



fc White| ( [19971 NFW) profile and (3) a local hot gas 
disk. The first profile is among the puffiest distributions 
one might consider, though it is in fact very similar to 
predictions for realistic galaxy-size halos with gaseous ha- 
los that are in both h ydrostatic and thermal equilibrium 
(Sharma et al. 2012). The second profile we consider, 
the MFW hot halo, is the most centrally concentrated 
distribution that one could envision for a hot gas around 
a galaxy, given that it is the profile that arises from the 
collapse of non-interacting dust. While it is unclear how 
such a profile could arise in a real astrophysical plasma, 
this is an assumption that often appears in the literature. 



1997 Bullock et al. 2001). The implied mass m miss 
ing baryons is approximately 10 11 M Q . We then explore 
the three aforementioned hot gas distributions, includ- 
ing two coronal models that contain Mhot = 10 11 M Q of 
material, in this context. 

2.1. Extended Adiabatic Halo: MB 

The first model we explore, MB, assumes that the hot 
halo is distributed as adiabatic gas with polytropic in- 
dex of 5/3 that is in hydrostatic equilibrium within the 
Milky Way's NFW dark matter halo of concentration C v 
( [Mailer fc Bulb^[2004| : 



pT(r)=Pv 



3.7 



Ln(l + x) 



3.7 



ln(l + C v ) 



3/2 



(!) 

where r is the radius from the Galactic center and x = 
r/R 8 , with R s the scale radius of the dark matter halo. 
The density p v is the gas density at the virial radius. 
In our fiducial case, we chose p v such that the integrated 
hot gas mass within R v is A/hot = 1O 11 M . The resulting 
temperature profile is 



/tiMB 
9 



(r) = T v \l + — ln(l + x) - ~ ln(l + C y ) 

" (2) 

where T v is the hot halo temperature at R v . We set T v by 
requiring that the X-ray emission-weighted temperature 
be consistent with XMM-Newton and Suzaku observa- 
tions, which both give halo temperatures consistent with 
Th = 0.2 keV, with errors of ~ 10%. The normalization, 
T v , is then determined by: 



T h = 



Jp2(r)A(T,Z g )T(r)dr 
Jpi(r)A(T,Z g )dr ' 



(3) 



where A(T , Z g ) is the cooling function (Sutherland & Do- 
pita|1993 ) and we assume a Z g = Q.3Zq (res ults tor dit- 
terent choices of Z g are explored in Section 4.1). The 
implied density, temperature, pressure, and mass pro- 
files for the MB model are shown as solid black lines in 
Figure [T] 

2.2. Cuspy Halo: NFW 

As mentioned in the introduction, it is common to 
make the simplifying assumption that the density dis- 
tribution of the hot gas around galaxies traces what is 
expected for the dark matter alone: 



Pa 

x(l + x) 2 ' 



(4) 



As before, x = r/R s , where R s = R v /C h is the NFW 
scale radius. We explicitly allow Cu to be different than 
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Fig. 1. — Hot gas profiles as a function of radius. Panel (a): density; panel (b): temperature; panel (c): pressure; panel (d): hot gas 
mass. The predictions from the MB model are plotted in solid black, the DISK model in dotted blue, and NFW models with C v = 12 (3) 
in red dash (dot-dash). For the DISK model, the x-axis is actually the vertical distance (z) from the disk. The square symbols in panel (a) 
show estimat es of the hot gas density req uired to explain the lack of HI in Milky Way dwarf galaxi es as derived from ram pressure s t riping 
arguments by |Grcevich fe Putman] ( |2009^ . The square symbols in panel (c) are derived by |Stanimirovic et al.||2002[ l and |Fox et al.||2005} 
under the assumption that high velocity clouds in the Milky Way halo are pressure-confined, bee text tor detausT 

C v of the background dark matter halo. As before, we 
fix the normalization p g by requiring that the total hot 
gas mass within the virial radius equals 10 11 M Q . For the 
NFW profile, we assume a constant hot halo temperature 
of Tfr. This requires an unusual equation of state profile 
that varies with radius for self-consistency, but we adopt 
it here for comparison to previous work. 

We explore two cases for the NFW hot halo concen- 
tration: C'h = 12 (which mimics the dark matter back- 
ground exactly: Ch = C v ) and also a low-concentration 
case with Ch = 3. The implied density, temperature, 
pressure, and mass profiles for the high (low) concentra- 
tion NFW case are shown by the dashed (dot-dashed) 
lines in Figure [T] 

2.3. Local Model: DISK 

For comparison, we also calculate the hot gas distri- 
bution in an exponential disk model, which h as been fa- 
vored by recent X-ray observations. Following Yao et al. 



T 9 DISK (z)=r exp 



z 

"hrfj ' 



(5) 



(6) 



Here £ is the volume filling factor and is assumed to be 
1, z is the vertical distance from the disk, po and To are 
the gas density and temperature at the disk mid-plane, 
and hp and hr are the scale height of the hot gas density 
and temperature distributions, respectively. 



1.4 x \Q~ A cm" 



2.3 kpc, T 



We set po 

10 6 4 K, and hr = 5.6 k pc based on observat i ons of the 
PKS2155-304 sight line j|Hagihara et al.|2010[ ) 
d2009ab observed the LMC X-3 sight 



Yao et al. 
1 reached 

similar conclusion. This model is plotted in blue in Fig- 
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ure[T] the y-axis is the vertical distance from the Galactic 
disk(|z|) since in the model the density and temperature 
depends only on this distance. 



3. OBSERVATIONAL CONSTRAINTS 

3.1. Indirect Probes 

One of the strongest pieces of indirect evidence for an 
extended hot gas reservoir around the Galaxy comes from 
the lack of gas detected in small dwarf satellite galaxies 
around the Milky Way. Dwarf galaxies tend to be gas rich 
unless they are within close (~ 300 kpc) proximity of a 
larger system, an observation that is usually interpreted 
as arising from ram-p ressure stripping |Lin fc Faber|1983 
IMoore fc Davisl[T994l 



Recently, |Grcevich fc Putman ( 2009 ) examined the HI 
content of the Local Group dwart galaxies. They found 
a cut-off radius of ~ 270 kpc around the Milky Way 
or Andromeda, below which no dwarf spheroidal galaxy 
contains detected HI gas. They argued the ram-pressure 
stripping produced by a hot, extended halo gas can ex- 
plain the lack of the HI gas in these dwarfs, and they 
constrained this gas density by the measurements of four 
dwarfs, Carina, Ursa Minor, Sculptor, and Fornax, be- 
tween ~ 20 to 100 kpc. The data points and their error 
bars are plotted in the top left panel of Figure m These 
measurements are much too extended to be explained by 
a local disk model (blue dotted) and too low density to 
be explained by the cuspy NFW models (red); they are 
matched reasonably well by the MB distribution, how- 
ever. 

A completely distinct indication for extended hot gas 
comes from studying gas clou ds around the galaxy . Rely 



ing on the Arecibo telescope, Stanimirovic et al. (20021 



observed the Magellanic Stream in Hi 21-cm emission 
They found that the most likely mechanism to confine 
the stream clouds is via pressure from a hot halo. They 
placed upper limits on the hot gas pressure at 10 3 and 
3 x 10~ 2 cm~ 3 K at 15 and 45 kpc, res pectively (square 



symbols in Figure n] bottom left panel). |Fox et al.|f2005| ) 
have reached similar conclusions when studying high ve- 
locity clouds (HVCs) in the halo. These authors con- 
cluded that if HVCs are pressure-confined, then the pres- 
sure of the surrounding medium must be approximately 
530, 140, and 50 cm~ 3 K, at distances of 10, 50, and 
100 kpc, respectively (triangles in Figure [TJ bottom left 
panel). Again, while neither the NFW nor DISK mod- 
els matches these constraints, the extended MB model 
provides pressure support at approximately the required 
level. 

Of course, neither the cloud-confinement arguments 
nor the ram-pressure stripping arguments individually 
demands the existence of an extended hot halo. It may 
be possible that some other environmental process can 
explain the lack of gas in local dwarfs, though there is no 
obvious candidate for such a mechanism. Similarly, it is 
possible that warm gas clouds in the Galactic halo are not 
pressure-confined. The confinement could, in principle, 
be gravitation al if the clouds are em bedded within dark 
matter halos ( [Sternberg et al.||2002 ). Nevertheless, it is 
interesting to recognize that an extended hot gas halo 
containing a substantial fraction of the halo's baryons is 
capable of matching both constraints quite well. 



3.2. X-ray Emission 

X-ray emission provides a strong constraint on the ex- 
tent and density of hot gas around the Galaxy. We 
consider data from two recent survey s with the XMM- 
Newto n and Suzaku X-ray telescopes. Henley fc Shelton] 
( |2010| studied 26 high-latitude XMM-Newton observa- 
tions of diffuse X-ray emission > 30° where b is the 
Galactic latitude). The Galactic longitudes (I) of these 
observations were selected to be between 120° and 240° 
to avoid enhanced emission from the Galactic center. 



Yoshino et al. (20091 studied the soft, diffuse X-ray emis- 
sion in 12 fields observed with the Suzaku X-ray Tele- 
scope. The Galactic longitudes and latitudes are similar 
to those XMM fields, i.e., high latitudes and away from 
the direction of the Galactic center, except two fields at 
b = 10° and 20°. 

While the exact details differ between the two analyses, 
the general procedures are similar (see their papers for 
details). After cleaning the spectrum, they fitted the 
observed spectrum to three components: (1) the cosmic 
X-ray background; (2) local thermal plasma emission at 
a temperature around 0.1 keV from solar wind charge 
exchange (SWCX) and local hot bubble (LHB) emission; 
and (3) distant thermal emission from near the Galaxy 
but beyond the LHB. The third component is the one 
that we consider here as likely arising from the distant 
halo. 

In the XMM data analysis, a single power law was 
adopted to model the cosmic X-ray background, so we 
use the results presented in Table 2 of Henley et al. For 
the Suzaku data, various models for the Cosmic X-ray 
background were tried and we present the the one that 
relied on a single power law (Yoshino et al.'s Table 4) 
in order to directly compare with the XMM result. We 
select the energy band between 0.4 and 2 keV (the default 
value adopted in the XMM data analysis) to compare the 
observed and modeled X-ray surface brightness (5° BS ). 
Since for the Suzaku data Sx BS were not calculated, we 
estimate the S£ BS using their models with the software 
package XSPEC; specifically, we use their model V (see 
their Table 4), which adopted a power law model of the 
CXB. The resultant normalized data is plotted in the left 
panel of Figure [2] The grey area indicates the full range 
of the data, which we will use as a benchmark comparison 
for our model predictions. 

We make predictions for the X-ray surface brightness 
from each of our hot gas models using 



(7) 



Here p e and pi are the electron and proton densities, and 
A(T) is the cooling function. In order to compare with 
observations, we use the Astrophysical Plasma Emission 
Code (APEC^to calculate the X-ray emissivity between 
0.4 and 2 keV as a function of plasma temperature. Since 
the X-ray emissivity is also a function of metal abun- 
dance, we assume a Z = O.SZ p for the MB and NFW 



models (Cen & Ostriker 2006 Rasmussen & Ponman 



2009 1 , and 1 Zp, tor the DISK modcl,~ where 'Zp, is the so 



Iar 

the effect of these choices on our predictions below 



bundance of Anders & Grevesse (1989). We explore 



5 See http://atomdb.org/ 
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Fig. 2. — X-ray emission constraints. Left panel: X-ray surface brightness between 0.4 and 2 keV as a function of Galactic latitude. 
XMM data are shown as squares, and Suzaku as diamonds (see text for details). The grey region shows the range of variation. Right panel: 
Model predictions of the X-ray surface brightness between 0.4 and 2 keV as a function of radius. Lines are the same as shown in Figure^ 
however, we plot the DISK model prediction for two b values, 0° (upper blue dotted curve) and 90° (lower blue dotted curve). The grey 
area is reproduced from the left panel. 



In the right panel of Figure [2j we plot the predicted 
surface brightness for our three models as the integrated 
emission coming from gas within a given distance (hor- 
izontal axis). The grey band is the constraint from ob- 
servations - the upper limits of the X-ray emission from 
the hot diffuse gas. Note that even though both XMM- 
Newton and Suzaku observations avoid the X-ray bright 
sky toward the Galactic center, we expect some sightlines 
may still experience local enhancement from supernova 
heating. 

Figure [2] clearly indicates the NFW profiles predict far 
too much X-ray emission. In particular, the C v = 12 case 
exceeds observations by more than two orders of magni- 
tude. This is co nsistent with the results of I Anderson fcl 



Bregman (20101, and it comes about because the density 
at the center 01 the NFW profile is far too high (see the 
left panel of Figure [l]). Indeed, the C v = 3 profile also 
substantially over-predicts the X-ray emission. We see 
from the rise of the line with distance that the majority 
of the X-ray emission comes from the inner ~ 10 kpc for 
both of the NFW models (red lines). Since the C h = 12 
hot gas NFW model is ruled out to an extreme degree, we 
will only include the Ch = 3 (low concentration) NFW 
halo in our comparisons for the remainder of this paper. 

In contrast to the NFW model, the MB profile (solid 
black) is fully consistent with the data, comfortably 
reaching the lower part of the observed region at ~ 
25 kpc. Recall that the MB profile has the same mass in 
hot baryons within R v as the NFW models, so the spa- 
tial distribution of the gas obviously plays a crucial role 
in determining the X-ray properties of the halo. The two 
dotted blue curves in Figure [2] represent the the DISK 
model predictions for in-plane (b — 0° ; upper curve) and 
vertical {b = 90°; lower curve) sight lines. The DISK 
model can also fit the observed value quite well at high 
Galactic latitude; however, at low Galactic latitude, it 
also predicts too much X-ray emission. This is partially 
a limitation of the simplified disk model assumed, which 
does not truncate in the radial plane. Unfortunately, 
with the current data we cannot tell whether there is a 



increasing trend at low b (see the right panel of Figure [l]); 
however, future targeted observations at low b may help 
distinguish between an extended hot halo like MB and a 
local hot distribution like the DISK model. 

3.3. Pulsar Dispersion Measure 

The dispersion measure (DM) of pulsars offers a direct 
probe of the electron distribution along the sight line 
towards a background pulsar at a distance D from the 
Earth: 



DM 



Pe(r)dr. 



(8) 



While most known pulsars live in the Milky Way disk, 
several have been discovered in the Large and Small Mag- 
ellanic Clouds (LMC and SMC), with a distance D ~ 50 
and 60 kpc, respectively. 
Two recent surve ys of pulsars in t he LM C a nd SMC 



were performed by Crawford et al. (2001) and Manch- 
ester et al. ( 2006 ) . They found a total of 21 pulsars, 
among which we selected 18 pulsars. Five of the selected 
pulsars are located in the SMC and the rest in the LMC. 
We discarded 3 pulsars with a DMsin|6| < 25 cm pc 
since the distribution of DMsin|6| of the Galactic pul- 
sars sugge sts those pulsars likely lie within the Milky 
Way disk ( |Crawford et al.||2001| |Gaensler et al.||2008[ ). 
We plot the DM of the remaining pulsars - which should 
be associated with the Magellanic Clouds - and their 
errors in Figure [3] For visual clarity, we randomly as- 
sign distances between 50 and 60 kpc for these pulsars. 
The pulsar DMs in LMC and SMC range from 40 to 
~ 200 cm -3 pc. Higher DMs for some of the pulsars 
indicate most likely a significant contribution from elec- 
trons in the LMC or SMC, so the lowest DM will provide 
an upper limit on the hot gas distributed between us and 
the LMC and SMC. 

There are two important potential sources for electrons 
giving rise to the DM. The first is the hot gas corona, 
which we will explore below. The second is the warm 
ionized medium (WIM) in and near the disk. It has long 
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Fig. 3. — Left panel: Pulsar dispersion measure as a function of distanc e. The dash-dotted lin e shows the contribution from the Warm 
Ionized Medium (WIM) in and near the disk using the parameterization of| Gaensler et al.| l |2008 ||. Otherwise the MB, NFW (C v = 3), and 
DISK line styles are the same as in Figure^ Only the corona distributed like an NP W provides a contribution larger than the WIM and 
is clearly ruled out by these data. Right panel: We also plot the total DM from the DISK and MB models (see the thick blue and dark 
curves, respectively). 

of ~100 kpc. The next brightest distant satellite (after 
the SMC and LMC) is the Fornax dwarf at a distance of 
~ 150 kpc. To our knowledge, there are no known pul- 
sars in Fornax ( McLaughlin fc Cordes|20"03 ), but further 
searches would be useful. 



been suggested that the distribution of the free electrons 
in th e WIM of the d isk follows a planar distribution (see, 
e.g., ( |Reynolds|p89| ): 



WIM 



„WIM 



exp 



H WIM 



(9) 



Here n^ IM and H™ m 



are the mid-plane WIM elec- 
tron density and the scale height, respectively. By fit- 
ting a total of 53 sightlines, Gaensler et al. (2008) found 



n™ M = 0.03li u: °u2 cm" 3 and tf„ wlM = 1010^ pc; we 
adopt these values here. The dash-dotted line in Fig- 
ure E] shows the contribution of this WIM distribution to 
the DM measure. We specifically calculate the contribu- 
tion along the line of sight toward the LMC. We see that 
it is significant, leaving very little room for additional 
electrons from a hot halo. 

The solid black, dashed red, and dotted blue lines in 
the left panel of Figure [3] show the DM coming from 
our MB, NFW (C h = 3)7 and DISK models, respec- 
tively. In each case we calculate the DM within a line- 
of-sight distance D pointed in the direction (I, b) via 
D 2 = r 2 + D 2 — 2rD a cos I cos b, where D = 8 kpc is the 
distance between the Sun and the Galactic Center and r 
is the Galactocentric distance. It is clear that both the 
MB and DISK models produce DM contributions that 
are well below the WIM, such that the total WIM+DM 
or WIM+DISK results are consistent with the data. The 
Ch = 3 NFW model, on the other hand, is clearly incon- 
sistent, especially when the total WIM+NFW contribu- 
tion is considered. 

Interestingly, the DM from the MB profile continues 
to rise to large radii. In principle, if we were confident 
about the contribution from the WIM, a dispersion mea- 
sure from more distant pulsars would provide interest- 
ing constraints such an extended profile. In the right 
panel of Figure |3]we show the detailed difference between 
the total DM ofthe MB+WIM model (black, solid line) 
and the DISK+WIM model (blue, dotted line). The two 
models already predict significant difference at a distance 



4. DISCUSSION 
4.1. Extended Halo: Parameter Dependence 

We have demonstrated that if one considers a hot halo 
profile with a low density (high entropy) core, as ex- 
pected for an adiabatic gas in hydrostatic equilibrium 
( Mailer &: Bullock"||2004[ ) or for a galaxy-size hot halo 
in hy drostatic and thermal equilibrium (Sharma et al. 
2012[ ) , then the Milky Way could in principle contain the 
universal baryon fraction within its virial radius. Specif- 
ically, in our fiducial MB model we considered a Milky 
Way halo of virial mass M v — 10 12 M and a hot halo 
of mass Mhot = 10 11 M , which would make the system 
baryonically closed. The X-ray surface brightness is quite 
sensitive to the chosen value of the metallicity of the gas; 
we have so far used Z g = 0.3 Zq. In this section, we in- 
vestigate how our MB halo results depend on our specific 
choice of virial mass, metallicity, and hot gas mass. 

The red dashed lines in the left panel of Figure HI show 
the density (top panel) and pressure (bottom panel) pro- 
files under the assumption that the Milky Way virial 
mass is at the upper end of the expected range, M v = 
2 x 10 12 M©, and that the system is baryonically closed 
(Mhot = 2.7 x 10 11 M Q ). Note that while the density 
itself appears to be consistent with the Grcevich & Put- 
man determination from HI gas stripping, the implied 
pressure profile is too high to explain the properties of 
gas clouds in the halo. The solid black lines show our 
fiducial model for reference and the blue, long-dashed 
lines show a similar halo that contains only ~ 2/3 of its 
baryons in total, with a hot halo mass of Mhot = 5x 10 10 
M©. This halo is fairly consistent with the Fox et al. 
pressure estimates, but it somewhat too low in density 
to explain the lack of HI in Milky Way dwarfs. 
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Fig. 4. — Left panel: Density (top) and pressure (bottom) distributions for MB profiles with Afhot = 2.7 (red dashed), 1.0 (solid black, 
fiducial), and 0.5 X 10 11 Mq (blue, long-dashed). Right panel: The X-ray surface brightness profiles for the same three models are shown 
with the same line types as in the left panel. In addition, the dotted and dot-dashed lines show our fiducial halo with both higher metallicity 
(Z g = IZq, top) and lower metallicity (O.IZq, bottom), respectively. 



The same group of models are presented in the right 
panel in comparison to the X-ray surface brightness con- 
straints for our fiducial metallicity assumption Z g = 0.3. 
They are all consistent with these data. The predicted 
Sx values are also sensitive to metallicity, as metals in- 
crease the cooling efficiency. To investigate the effect 
of metals, we calculate Sx for two additional mctallici- 
ties within our fiducial MB halo: Z g = 1 Zq (the dotted 
dark line in Figure[4]), and Z g = 0.1 Zq (the dashed dark 
line). Although the predicted Sx is still consistent with 
observations for Z g = 1 Zq , it is on the high end of the 
observed Sx- 

5. SUMMARY 

Whether or not galaxies like the Milky Way host sub- 
stantial masses of baryons in hot coronae has a funda- 
mental impact on our understanding of galaxy formation 
and evolution. In this paper, we examined the latest X- 
ray emission and pulsar dispersion measure data to test 
for the presence of hot baryons in the distant halo of our 
Milky Way. The three models that we used to describe 
the hot gas distribution in the halo are: (1) a NFW-type; 
(2) an extended hot halo motivated by MB04; and (3) a 
disk distribution. 

We found that for a baryonically closed Milky Way, 
the hot gas cannot follow an NFW profile, either for 
a standard halo concentration (Ct — 12) or for a low 
concentration {Ch = 3): the NFW profile predicts both 
too much X-ray emission and a pulsar dispersion mea- 
sure that exceeds the latest observational data. The 
baryon fraction in such models must be substantially 
low er (fa ~ 0.01 — 0.02) to be c onsistent with data (see 



also Anderson & Bregman 



2010) 



The other two classes of models we have considered 
- extended halos of hot gas in hydrostatic equilibrium 
(Mailer & Bullock 2004) and a hot disk of gas - can 
both be made consistent with existing X-ray emission 
and pulsar dispersion measure data. These two models 
predict very different properties for the hot gas content 
of the Milky Way, however: in the DISK model, the hot 
gas contributes a negligible amount (~ 10 8 M Q ) to the 
Milky Way's baryon budget, while the MB extended halo 



model can contain enough hot gas to make the Milky 
Way baryonically closed. 

As we were completing this work, Gupta et al. (2012) 
presented a complementary study of the hot gas con- 
tent of the Milky Way halo using X-ray spectra of back- 
ground AGNs. Based on the detected z = absorption 
lines produced by highly ionized oxygen, O VII and O 
VIII, and a joint analysis with the Galactic halo emission, 
Gupta et al. argued for the existence of an extended hot 
gas around our Galaxy, with a radius of over 100 kpc 
and a total mass in excess of 10 10 Mq. Their conclusion 
is fully consistent with our results for extended hot gas 
profiles, and incorporating X-ray absorption data, along 
with considerations of uncertainties associated with the 
gas metallicity and ionization mechanism(s), would likely 
be a fruitful avenue for future extensions of our current 
work. 

Incorporating additional indirect constraints - i.e., the 
lack of gas in all Milky Way dwarf spheroidals and the 
possibility that high velocity clouds are pressure-confined 
by a hot ambient medium - may indicate that a MB-type 
model with an extended hot halo is favored. This inter- 
pretation should be testable in the near future, as obser- 
vations of X-ray emission at low latitudes < 20°) 
can distinguish between the MB and the DISK mod- 
els. Such observations ideally should be performed at 
90 < I < 270° to avoid the contamination from the 
Galactic center. If pulsars can be detected in distant 
dwarf spheroidal galaxies with recent star formation - 
e.g., Fornax and Leo I - dispersion measures for these 
objects would also be very constraining for the distribu- 
tion of hot gas around the Milky Way. 
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